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c^ ; ABSTRACT 

Vh ' Effects of Coulomb interactions and lattice fluctuations in the optical ab- 

^ • sorption spectra of Ceo and C70 are theoretically investigated by using a 

C . tight binding model with long-range Coulomb interaction and bond disorder. 

Qj ! Anisotropy effects in C70 are also considered. Optical spectra are calculated 

by using the Hartree-Fock approximation followed by the configuration in- 
k^ teraction method. The main conclusions are as follows: (1) The broad peaks 

'j_| ■ at excitation energies, 3.7eV, 4.7eV, and 5.7eV, observed in experiments of 

"^ Ceo molecules in a solution are reasonably described by the present theory. 

Peak positions and relative oscillator strengths are in overall agreement with 
the experiments. The broadening of peaks by lattice fluctuations is well sim- 
ulated by the bond disorder model. (2) The optical gap of C70 is larger when 
the electric field of light is parallel to the long axis of the molecule. The shape 
of the frequency dispersion also depends on the orientation of the molecule. 
These properties are common in the free electron model and the model with 
Coulomb interactions. (3) The spectrum of C70 averaged over bond disorder 
and random orientations is compared with experiments in a solution. There 
is an overall agreement about the spectral shape. Differences in the spectra 
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of Ceo and C70 are discussed in connection with the symmetry reduction from 
a soccerball to a rugbybaU. 

PACS numbers: 78.66. Qn, 78.20.Dj, 71.35.+Z, 31.20.Tz 
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I. INTRODUCTION 

Recently, fullerenes Cat (A^ = 60, 70, 76, and so on) with hollow cage struc- 
tures have been intensively investigated. Many optical experiments have 
been performed, and interesting properties originating from vr electrons de- 
localized on molecule surfaces have been revealed. They include the optical 
absorption spectra of Ceo [1;2] and C70 [1], and the large optical nonlinearity 
of Cqo [3,4]. The nonlinearity in the third harmonic generation (THG) is of 
the order lO^^^esu, and this largeness is attractive in the scientific as well 
as technological interests. The C70 thin films show a large THG magnitude 
extending to 10~^^ esu also [5]. Therefore, the fullerene family is generally 
attractive for possible application to nonlinear optical devices in near future. 

In order to analyze the optical properties and to clarify mechanisms of 
the large nonlinearity, we have studied the linear absorption and the THG 
of Ceo by using a tight binding model [6] and a model with a long-range 
Coulomb interaction [7]. A free electron model yields the THG magnitudes 
which are in agreement with the experiment of Ceo [2], while the calculated 
linear absorption spectrum is not in satisfactory agreement with experiments 
[6]. If Coulomb interactions are taken into account, the absorption spectra 
are in overall agreement with the experiment, although the magnitude of the 
THG decreases by a factor about 0.1 [7]. 

The first purpose of this paper is to report the details of the analysis of 
optical absorption spectra of Ceo- In the previous brief paper [7], the effects 
of the long ranged Coulomb interactions have been considered in the tight 
binding model for a single molecule. The model has been analyzed by the 
restricted Hartree-Fock approximation followed by the single excitation con- 
figuration interaction method (single CI) [8] . We use the same method in this 
paper. In the experiments of actual materials, the optical spectra become 
broad mainly due to the lattice fiuctuations. We are going to simulate the 
effects by using the bond disorder model with Gaussian distribution. We will 
take a sample average over disorder potentials. This method has been used 
in order to describe the optical properties of conjugated polymers [9,10]. We 
will show that main features observed in experiments of Ceo molecules in 
a solution [2] are reasonably described by the present theory. Speaking in 
detail, peak positions in the frequency dispersion of the spectra agree well 
with the experiment. The relative oscillator strengths are in overall agree- 
ment with the experimental data. Furthermore, the broadening of absorption 
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peaks by lattice fluctuations is well simulated by the bond disorder model 
with a reasonable disorder strength. 

Secondly, we extend the calculations to one of higher fullerenes: C70. The 
molecular shape becomes elongated in going from Cgo to C70. Then, we ex- 
pect that the optical response varies depending on the molecular orientation 
relative to the electric field of light. We will discuss the anisotropy for the free 
electron model and for the model with the Coulomb potentials. It is found 
that the optical gap is larger when the electric field is parallel to the long axis 
of the C70 molecule. The line shape of the absorption spectrum is anisotropic, 
too. These properties do not depend on whether the Coulomb interactions 
are present or not. We expect that the anisotropy will be observed in the 
solids of C70. 

Finally, we consider a situation in which the C70 molecule is rotated 
randomly and the anisotropy of the optical absorption is averaged out. The 
optical spectrum calculated by the model with Coulomb interactions and 
disorder potentials are compared with the experimental data [1] for molecules 
in a solution. We will conclude that there is an overall agreement. The 
absorption around 2.7 eV in the experiment can be interpreted as an allowed 
transition. There are several allowed transitions below the energy 2.7eV, but 
their oscillator strengths are smaller. We can relate the spectra of Ceo and 
C70 by considering the symmetry reduction from a soccerball to a rugbyball. 

We shall report calculations in the following way. In the next section, our 
model and its meaning are explained. The calculation method of the single 
CI is also summarized. In Sec. Ill, we report about Cgo and compare with 
the experiment. In Sec. IV, we describe optical spectra of C70. The paper is 
closed with several remarks in the final section. 



II. MODEL AND FORMALISM 

In order to consider optical spectra of Ceo and C70, we use the following 
hamiltonian: 

H = Hq + H],ond + Hmt. (2.1) 



The first term Hq of Eq. (2.1) is tlie tiglit binding model. We use tlie 
following form for Ceo^ 

Ho= Y: °(-tD)(4c,,. + h.c.) + Y: '(-ts)(4c„. + h.c.), (2.2) 

where Cj^o- is an annihilation operator of a vr-electron at site i with spin a; the 
sum with the symbol D (or S) is taken over all the pairs (i, j) of neighboring 
atoms with a double (or single) bond; and t-o = t+(2/3)t' and ts = t — {l/3)t' 
are the hopping integrals. The mean hopping is taken as t. In this paper, we 
use the difference between tj) and t^, t' = O.lt, as a typical value. The choice 
of t' does not affect on the results so strongly, because main contributions 
come from the strong Coulomb potential. We use the next form for C70: 

Ho = -tT. clc,,^, (2.3) 

assuming a constant hopping integral. We could use a realistic combination 
of hoppings, but the results do not change so much like in the Cgo case. 
Effects of zero point vibrations and thermal fluctuation of the lattice are 
modeled by the bond disorder model which is the second term of Eq. (2.1): 

-f^bond = Yl (^hj4,aCj,a- (2.4) 

Here, 6tij is the disorder potential at the bond {i,j). We can estimate 
the strength of the disorder (standard deviation) tg from the results [11] by 
the extended Su-Schrieffer-Heeger (SSH) model [12]. The value would be 
ts ~ 0.05 — O.lt. This is of the similar magnitude as in the fullerene tubules 
and conjugated polymers [13]. We shall treat interactions among vr-electrons 
by the following model: 



Y W^Ki)(IIcl,^Ci,^ - l)(5]c;-^c,- 



J 



where rij is the distance between the ith and jth sites and 

W(r) = ^ (2.6) 

y'(l/f/)2 + {r/roVy 



is the Ohno potential. The quantity U is the strength of the onsite interac- 
tion, V means the strength of the long range Coulomb interaction, and tq is 
the average bond length. 

The model is treated by the Hartree-Fock approximation and the single 
CI [8]. After the Hartree-Fock approximation H =^ H-^p, we divide the total 
hamiltonian as iJ = i^HF + H'. The term H' becomes 

H' = f^E(4T^^,T-P.T)(4i^U-Pu) (2.7) 

i 
+ Y. ^(^ij)E(4,aCi,a - Pi,a){c]^^Cj,r - Pj,r) 

"T~ / A ^i,i,cr^i.(T^i,<T ' 'Tj,i,crCi^(jCj^cr Ti,j,aTj,i,a)\i 
a 

where pi^a = (cI,o-Cj,cr) and tij^o- = (cI,o-Cj>)- When we write the Hartree-Fock 
ground state \g) = IlAioccupicd Ca,|CJ^,JO) and the single electron hole excita- 
tions l/iA) = c|j^^CA,r|fi') {p means an unoccupied state; we assume both singlet 
and triplet excitations in this abbreviated notation), the matrix elements of 
the Hartree-Fock part and the excitation hamiltonian become as follows: 

{p'X'\{HnF-{H^F))\pX) = V,A',a(^m-^a), (2.8) 

{p'\'\iH'-{H'))\p\) = 2J6s-K, (2.9) 

where E^j, is the energy of the Hartree-Fock state, Ss = I for spin singlet, 
Ss = for spin triplet, and 



J(/i', A'; /i. A) = E^M(p'K)(A'K)(j|/i)(j|A) (2.10) 

Kip',\';p,\) = ^\/,,(^'|j)(A'K)(j|/i)(^|A) (2.11) 



*j 



with Vi^i = U, Vij = W{rij) for i ^ j. The diagonalization of the total 
hamiltonian H gives the set of the excited states {\k,)} within the single CI 
method. In the actual calculation, we limit the spin configurations to the 
singlet excitations. 

The absorption spectrum in the case of the electric field parallel to the x 
axis is proportional to 

J2E,PiLU-E,){g\x\K){K\x\g). (2.12) 



Here, P{uj) = ^/[n^u^ + 7^)] is the Lorentzian distribution (7 is the width), 
Ei^ is the electron-hole excitation energy, and \g) means the ground state. 

In this paper, all the quantities with the energy dimension are shown 
in the units of t. We varied the parameters of Coulomb interaction within 
< V < U < 5t, and we report here the two representative cases: U = V = 
and U = 2V = 4t. Results in the latter case turn out to be in overall 
agreement with experiments. 



III. OPTICAL ABSORPTION IN C 
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This section is devoted to optical spectra of Cgo, the structure of which is 
shown in Fig. 1(a). The calculated spectra are nearly isotropic independent 
of the direction of the electric field against the molecular orientation. 

A. Coulomb interaction effects 

The dispersion of the optical absorption is calculated for the free electron 
case U = V = 0, and displayed in Fig. 2(a). The Lorentzian broadening 
7 = 0.0 It is used. The distribution of the relative oscillator strength and the 
positions of the peaks agree with those in the previous calculations [6]. The 
second peak at about 1.2t has the largest oscillator strength. 

Fig. 2(b) shows the absorption for U = 4t and V = 2t. There are several 
peaks, among which the peak at ~ 2.6t has the largest oscillator strength. 
This is in contrast to the free electron case {U = V = 0), where the oscillator 
strength is the largest at around 1.2t. The Coulomb interaction tends to 
reduce the oscillator strengths of the lower-energy peaks and to enhance 
those of the higher-energy peaks. In addition, it shifts overall peak positions 
to higher energies. Qualitatively similar results have been reported in terms 
of the random phase approximation [14]. The overall distribution of the 
oscillator strengths in Fig. 2(b) is in good agreement with experiments [2] 
when we assume t = 1.8eV. The origin of the large broadening of the peaks 
may be lattice fluctuations. This is to be discussed in the next subsection. 

Now, we can estimate a dielectric constant e. The magnitude of the 
Coulomb interaction at the mean carbon bond length is W^(ro) = UV/\^U^ + V^. 
This gives an estimate: e = e^/[roW^(ro)] = 3.11. The long distance limit 
W{r ^ ro) ~ Vro/r gives e = 2.80. The independent study of the lifetime 



of positrons in Ceo gives the high frequency dielectric constant ^oo ~ 3.5 [15]. 
Therefore, the present estimation of the dielectric constant near the value 3 
would be reasonable. 

B. Lattice fluctuation effects 

Here, we should like to take into account of the lattice fluctuation effects. We 
simulate the effects by the bond disorder model with a Gaussian distribution. 
The application of the model has been successful in describing the lattice 
fluctuation effects in conjugated polymers [9]. The estimation by the results 
of the extended SSH model for Ceo gives the order tg ~ 0.05 — O.lt [11]. 
We change the disorder strength around these values. Averaging over 100 
samples turned out to be enough to obtain the smooth optical spectra. 

Fig. 2(c) shows the calculated optical spectrum for tg = 0.09t. For 
comparison we show in Fig. 2(d) the experimental data obtained by Ren et 
al. for Ceo in a solution and a Ceo thin film [2]. The data are shown for 
comparison assuming t = 1.8eV. Due to the broadening, the peaks near the 
energy 1.8t in Fig. 2(b) merge into a single broad feature in Fig. 2(c). The 
several peaks around 3.4t obtain the same effect as well. Overally, there are 
three main features in the energy dependence of the calculated data. Their 
positions in the energy agree well with the experimental data, 2.0t, 2.6t, and 
3. It, obtained from the solution. The relative oscillator strengths fairly agree 
also. Furthermore, the widths of the features can be well simulated by the 
broadening from the bond disorder model. Experimentally, the peaks would 
be broadened by various kinds of phonon modes: intermolecular libration and 
intramolecular vibrations. They give rise to many sidebands, which broaden 
into several main features. The disorder strength, adopted here, simulates 
well the total broadening as the net contribution from the various phonon 
modes. 

In the solid, the broadening is larger and a broad hump appears around 
the energy 2.8eV (= 1.5t), as is shown in Fig. 2(d). There is a tiny structure 
at 3eV in the corresponding region of the absorption spectrum of Ceo in a 
solution (Fig. 2(d) and Ref. 17). We assume that the forbidden transitions 
in that energy region become partially allowed due to lattice fluctuations or 
intermolecular interactions. If the lattice fluctuations are effective, the effect 
can be simulated by the bond disorder model. In Fig. 2(c), the absorption 
in the low energy part multiplied by the factor 10 is also shown. The several 



forbidden transitions around 1.2t — 1.6t become allowed by the disorder, 
giving rise to an absorption tail in the lower energy region. The relative 
magnitude is similar to the solution data shown in Fig. 2(d). This might 
be also the origin of the 2.8eV hump in the solid. However, the strength 
of the absorption relative to the main peaks in our theory is about one 
order of magnitude smaller than that in the experiment. If we increase 
disorder drastically to explain the magnitude, we get absorption peaks that 
are too broad compared with the experiments. This is true for all the cases 
of Coulomb interaction parameters we tried. Therefore, it seems difficult to 
explain the oscillator strength of the partially allowed transitions in solids 
quantitatively. Additional effects such as intermolecular interactions will be 
necessary, and this problem will be investigated in a separate paper [16]. We 
note that the Cgo oxide shows a weak absorption peak in a similar energy 
region [18]. This could also be interpreted by using a model with impurities. 

IV. OPTICAL ABSORPTION IN C70 

A. Anisotropy effects 

The C70 molecule has a rugbyball-like shape. Figs. 1(b) and (c) display the 
molecule. We locate the long axis of the molecule along the x direction. Fig. 
1(b) shows the molecule which is projected along the z axis. The x axis goes 
from the left to the right of the figure. Fig. 1(c) is the molecule projected 
along the x axis. We expect that the absorption becomes anisotropic with 
respect to the direction of the electric field of light. We look at anisotropy 
effects in the free electron model as well as in the model with Coulomb 
interactions. 

Figs. 3 (a), (b), and (c) show absorption spectra of the free electron case 
when the electric field is in the x, y, and z directions, respectively. Figs. 3 
(b) and (c) have the almost same spectral shape. But, the spectral shape is 
very different when the electric field is parallel to the long axis. The optical 
gap is about 0.7t in Fig. 3(a), while it is about 0.6t in Fig. 3(b) and (c). 
This accords with the result of the recent Hiickel calculation [19] in that the 
optical gap of C70 is larger when the electric field is parallel to the long axis 
of the molecule. 

Next, we shall turn to Coulomb interaction effects. We display the optical 
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spectra for U = 4t and V = 2t. We will later find that these parameters are 
appropriate for C70 also. Figs. 4 (a), (b), and (c) show the spectra for the 
respective molecular orientations corresponding to Figs. 3(a-c). It seems 
that there is a strong anisotropy when the Coulomb interactions are present, 
too. The optical gap in Fig. 4(a) is about l.lt with a weak absorption peak. 
It is about l.Ot in Figs. 4 (b) and (c). Therefore, the optical gap becomes 
larger when the electric field is parallel to the long axis of the molecule. This 
property does not change whether there are Coulomb interactions or not. We 
highly expect that the anisotropy will be clearly observed in the solid of C70. 

B. Orientationally Averaged Spectra 

1. Coulomb interaction effects 

The optical absorption data have been reported for molecules in solution [1]. 
In order to compare the present theory with the experiment, we have ro- 
tated the molecule randomly and averaged the spectra over 100 orientations. 
Averaged data without bond disorder are shown in Figs. 5(a) and (b). We 
use the Lorentzian broadening 7 = O.Olt as before. Fig. 5(a) is for the free 
electron case, and Fig. 5(b) is for the Coulomb interaction strengths U = 4t 
and V = 2t. We can easily see that the peaks in three figures of Fig. 3 
remain in Fig. 5(a). The same thing can be said in the relation between Fig. 
4 and Fig. 5(b). The effects of the Coulomb interactions are similar to those 
of the Cgo case: the Coulomb potentials tend to shift peaks to higher ener- 
gies, and oscillator strengths at higher energies have relatively larger weight 
than in the free electron case. The shift is due to the fact that the inter site 
Coulomb interactions give rise to bond order parameters in the Hartree-Fock 
approximation. The bond order effectively enhances the intersite hopping 
integrals, and thus the energy gap. 

2. Lattice fluctuation effects 

Now, we treat lattice fluctuation effects by using the bond disorder model. 
As we expect that the width of lattice fluctuations in C70 is similar to that 
in Ceo even if the molecular structure is different, we assume the same tg = 
0.09t in this calculation again. The orientation of the molecule and disorder 
potentials are changed 100 times, and the optical spectra are averaged. 
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Fig. 5(c) shows the calculated optical absorption. Fig. 5(d) shows the ex- 
perimental data of C70 in a solution taken from [1]. The abscissa is scaled by 
using t = 1.9eV. The optical gap near l.Ot accords well with the experiment. 
The small peaks around 1.8t may correspond to those of the experiment in 
the similar energy region. Commonly, there is a dip of the optical absorption 
at about 2.0t in the calculation and experiment. There is a maximum at 2.8t 
in the two data. Therefore, we conclude that there is an overall agreement 
between the calculation and the experiment. The quantum lattice fluctuation 
effects are also simulated well by the bond disorder model. 

V. CONCLUDING REMARKS 

The most striking difference between optical spectra of Cqq in solution and of 
Ceo solids is that the feature around 2.8eV has the strong oscillator strength 
only in solids. In the present calculation for U = 4t and V = 2t without 
disorder, there are dipole forbidden transitions at the corresponding energy. 
The disorder makes these transitions partially allowed. We have discussed 
that this effect is one of the origins which make the forbidden transitions 
observed. However, the theoretical oscillator strengths seem to be smaller 
by an order of magnitude than that in the experiments. We need further 
study in order to solve this problem. For example, we are now studying 
intermolecular interaction effects in addition to lattice fluctuations. This 
will be reported separately [16]. 

The absorption spectra of C70 become highly anisotropic depending on the 
orientation of the molecule relative to the polarization of light. Now, single 
crystals of C70 have been made and their optical data have been reported 
[5]. Therefore, an extension of the present calculations to C70 solids is quite 
interesting. We expect that the anisotropy of the optical spectra of the 
molecule will remain in solids, too. 

In view of the agreement of our theory with experiments in Cqq and C70, 
it is very interesting to look at optical properties and excitonic effects of 
higher fullerenes where the number of carbons is larger than 70. The method 
of the single CI should be valid enough in these systems, too. The number 
of isomers rapidly increases as the carbon number increases as 76, 78, 82, 
84, and so on. The calculated optical data will reflect the electronic and 
lattice structures of isomers, and they will be mutually different. Thus, the 
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calculations will be useful in order to distinguish isomers by using optical 
experiments. 
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FIGURE CAPTIONS 

FIG. 1. Molecular structures for (a) Cqo, and (b,c) C70. In (b), the long 
axis {x axis) penetrates from the left to the right of the figure. The mirror 
symmetry planes are perpendicular to the figure. In (c), the C70 is projected 
along the long axis. 

FIG. 2. Optical absorption spectra for Ceo shown in arbitrary units. The 
abscissa is scaled by t. The spectra are calculated with the parameters (a) 
U = V = 0, (b,c) U = At and V = 2t. The Lorentzian broadening 7 = 
O.Olt is used. In (c), lattice fluctuations are taken into account by the bond 
disorder of the strength tg = 0.09t. (d) The experimental spectra (Ref. [2]) 
of molecules in a solution (solid line) and of a Cgo film (dotted line). We use 
t = 1.8eV. 

FIG. 3. Optical absorption spectra for C70 shown in arbitrary units. The 
abscissa is scaled by t. The interaction strengths are U = V = 0. The 
Lorentzian broadening 7 = O.Olt is used. In (a), (b), and (c), the electric 
field is parallel to the x, y, and z axes, respectively. 

FIG. 4. Optical absorption spectra for C70 shown in arbitrary units. The 
abscissa is scaled by t. The interaction strengths are U = 4t and V = 2t. 
The Lorentzian broadening 7 = O.Olt is used. In (a), (b), and (c), the electric 
field is parallel to the x, y, and z axes, respectively. 

FIG. 5. Orientationally averaged optical absorption spectra for C70 shown 
in arbitrary units. The abscissa is scaled by t. The spectra are calculated 
with the parameters (a) U = V = 0, (b,c) f/ = 4t and V = 2t. The molecule 
is rotated randomly in order to average the anisotropy. The Lorentzian 
broadening 7 = O.Olt is used. In (c), lattice fiuctuations are taken into 
account by the bond disorder of the strength tg = 0.09t. (d) The experimental 
spectrum (Ref. [1]) of molecules in a solution. We use t = 1.9eV. 
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